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INTRODUCTION: 

It  is  well  known  that  oxygen  has  significant  impact  on  cancer  treatment.  Our  hypothesis  for  this  research 
was:  (a)  Tumor  blood  vessels  were  leaky  and  therefore  acoustic  vibration  can  be  used  to  modulate  the 
leaky  vessels  and  induce  oxygenation  changes;  (b)  The  oxygenation  changes  can  be  detected  by  optical 
measurements.  Preliminary  studies  with  5  tumor-bearing  rats  demonstrate  that  ultrasonic  vibrations  can 
eidier  generate  significant  effects  (early  stage  tumors)  on  optical  measurements  or  no  effects  on  optical 
measurements  (late  stage  tumors).  During  the  past  year,  the  PI  has  devoted  his  efforts  on  quantify  both 
acoustic  vibration  and  optical  measurement  of  oxygenation. 

BODY: 

1.  Optimization  of  ultrasound  system  for  quantification  of  mechanical  stress  generated  by  the 
ultrasound  vibration  (in  task  2) 

The  ultrasound  system  that  the  PI  used  to  conduct  the  preliminary  animal  experiments  was  transmission 
only.  The  acoustic  pulses  of  various  energy  levels  were  delivered  to  the  tumors  and  the  acoustic  induced 
effects  were  measured  optically.  However,  the  optical  measurements  were  not  directly  proportional  to 
the  acoustic  energy  and  were  significantly  modulated  by  the  tumor  elastic  condition,  liiis  is  probably 
why  we  observed  significant  changes  for  two  early-stage  tumor  rats  but  no  changes  for  two  late  stage 
tumor  rats. 

To  directly  measure  ultrasonic  pressure  effects,  we  have  constructed  the  64-channel  receiving  ultrasoxmd 
system.  The  system  is  designed  for  real-time  RF  data  acquisition,  which  can  be  used  to  calculate  the 
mechanical  strain  generated  at  tumor  spot. 

Analog  Analog  Analog  Analog 


Fig.  1.  Block  diagram 
of  a  single  PCI  DAQ 
card 


The  64-channel  data  acquisition  system  includes  two  identical  standard  full-length  PCI  DAQ  cards  for 
personal  computers.  Each  PCI  DAQ  accepts  32-channel  analog  inputs  and  digitizes  all  the  inputs 
simultaneously  at  a  rate  of  40MSPS,  lObit.  The  block  diagram  of  a  single  PCI  DAQ  card  is  shown  in 
Fig.  1.  For  each  PCI  DAQ  card,  32-channel  analog  inputs  are  digitized  by  16  dual-channel  40MSPS, 
lObit  ADCs  MAXI  186  (Maxim  Inc).  Digitized  RF  data  are  stored  in  a  160bit  wide,  lOOMHz  FIFO 
memory  array.  The  RF  data  can  be  transferred  to  host  computer  at  a  data  transfer  rate  of  132MB/s 
through  a  PLX9054  (PLXTECH,  Inc)  PCI  interface  chip  when  operating  at  DMA  mode,  or  to  external 
DSP  board  at  a  data  transfer  rate  of 400MB/s  for  process.  Three  Altera  FPGAs  are  used  to  interface  with 
PCI  interface  chip,  external  DSP  board,  external  digital  I/O  lines  and  external  clock  and 
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synchronization.  An  LVDS  (Low  Voltage  Differential  Signaling)  clock  generator  and  distributor 
network  are  implemented  to  provide  very  accurate,  noise-immime  and  jitter-free  clocks  to  all  ADC 
channels.  In  PCB  (printed  circuit  board)  layout,  cares  are  taken  to  provide  a  constant  impedance,  equal- 
length  traces  of  ADC  clock  signals  to  ensure  accurate  phase  uniformity  when  digitizing  all  32  analog 
channels.  Proper  analog/digital  ground  separation  and  wide  band  input  RF  transformers  are  used  to 
provide  groxmd  isolation  and  reduce  digital  noise  and  cross-talk  between  analog  channels.  The 
photography  of  the  data  acquisition  card  is  shown  in  Fig.2.  Currently,  the  system  hardware 

construction  and  testing  are  completed  and  the  efforts  are  devoted  to  software  debugging  and 
implementing  strain  estimation  algorithm. 

Fig.2.  Photo  of  the  data 
acquisition  card  that  will  be 
used  to  estimate 
ultrasonically  generated 
stain,  which  may  induce 
tumor  oxygenation  changes 


2.  Optimization  of  optical  wavelength  selection  for  in  vivo  oxygenation  estimation  (in  tasks  1  and 
3) 

The  optical  measurements  observed  in  preliminary  rat  tumor  studies  were  simply  amplitude  changes  and 
quantitative  oxygenation  calculation  is  needed  to  evaluate  the  proposed  methods. 

Literature  data  regarding  to  optimal  optical  wavelength  for  tumor  hemoglobin  concentration  and  oxygen 
calculation  are  very  confusing.  Many  research  groups  have  constructed  multi-wavelength  systems  using 
either  frequency  domain  or  time-resolved  approach  [1-2].  The  total  number  of  wavelengths  deployed  in 
the  reported  systems  has  been  increased  dramatically  in  rec^t  years  varying  from  3  to  9  optical 
wavelengths.  With  more  wavelengths,  the  wavelength-dependent  tumor  absorption  changes  can  be 
obtained.  However,  the  computation  of  total  hemoglobin  and  hypoxia  distributions  requires  an  optimal 
set  of  wavelengths.  Fiuthermore,  the  data  acquisition  time  will  increase  dramatically  with  the  increasing 
wavelengths.  Recently,  the  PI  worked  with  other  Ph.D.  students  in  the  Lab  have  constructed  a 
frequency  domain  system  consisting  three  wavelengths  of  660nm,  780  nm  and  830  nm  [3]. 

Are  these  three  wavelengths  good  enough  to  accurately  and  robustly  compute  the  hemoglobin 
concentration  and  oxygenation  from  optical  measurements?  The  PI  has  analyzed  this  problem  from 
both  theory  and  practice.  The  analysis  is  as  follows. 


By  assuming  that  the  major  chromophores  are  deoxygenated  (deoxyHb)  and  oxygenated  (oxyHb) 
hemoglobin  in  the  wavelength  range  studied,  we  can  estimate  deoxyHb  and  oxyHb  concentrations  at 
each  imaging  voxel  by  inverting  the  following  equations 


p\  ph 

V 

deoxyHb{r') 

p^  p^ 

A 

oxyHb{r') 
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and  obtain  deoxyHb  and  oxyHb  as: 
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oxyHb{r') 


1 

— p^ 

^Hb02 

X 

4(''’) 

A 

1 

.4  ('■’). 

(2) 


The  total  hemoglobin  concentration  totalHb{r')-  deoxyHbir')  +  oxyHb(r')  and  oxygenation  saturation 

Y%  = - ^ - 100%  can  be  calculated  as: 

oxyHb{r' )  +  deoxy{r' ) 


and 
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where  ^  =  £^£^02  ^  principle,  any  two  wavelengths  in  the  NIR  window  can  be  used  to 

compute  the  total  hemoglobin  concentration  from  equation  (3).  However,  positive  coefficients  of 
and  in  equation  (3)  are  critical  to  guarantee  that  the  estimated  total  hemoglobin 

concentration  at  each  voxel  is  positive,  especially  in  the  background  region  where  the  absorption 
coefficients  are  small  and  consequently  the  relative  estimation  errors  of  absorption  coefficients  are 
larger  that  those  in  the  lesion  region.  The  second  practical  consideration  is  the  sensitivity  of  the  total 
hemoglobin  concentration  to  estimation  errors  of  (r')  and  The  sensitivity  is  defined  as  the 

partial  derivatives  of  totalHb  to  (r' )  and  (r* ) ,  which  are  given  as: 


dtotalHbir')  1  r  , 


DtotalHbir') 


(5) 


If  a  wavelength  after  isobestic  point  is  chosen  as  A^,  for  example  830  nm,  the  coefficient 
“^4)}  of  /^^(^')  is  positive  as  shown  in  Fig.  3  (circled  line),  however,  the  coefficient 

■^{{^4.  of  (^')is  positive  only  for  wavelength  /i,  <798  nm  (  *  line).  The  extinction 

coefficients  used  for  calculation  were  obtained  from  Ref.  [ph.d  dissertation]  in  a  natural  logarithmic 
scale  with  units  of  inverse  millimoles  times  inverse  centimeters.  However,  the  sensitivity  of  total 
hemoglobin  estimation  to  for  ^  <  780  nm  is  much  higher  than  that  of  (r' ) .  The  best  equal 

sensitivity  wavelength  range  of  for  given  =  830  nm  is  between  780  nm  and  785  nm. 
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Similarly,  any  two  wavelengths  in  the  NIR  window  can  be  used  to  compute  the  oxygenation  saturation 

from  equation  (4).  However,  the  sensitivity  of  Y%  to  ratio  is  significantly  different  for  different 

wavelength  pairs.  Figure  4  plots  the  Y%  vs.  the  ratio  for  wavelength  pairs  (660nm,  830nm,  circles), 
(710  nm,  830nm,  starts)  and  (780nm,  830nm,  crosses).  For  the  first  two  pairs,  the  Y%  lies  in  the 

— - >  0.2.  However,  for  the  wavelength  pair  of  (780  nm,  830nm),  the  Y%  is 


physiology  range  for 


greater  than  100%  for  - <0.7  and  reduces  quickly  to  zero  for  ratio  >  1.4.  Therefore,  any 

estimation  error  on  can  cause  the  Y%  out  of  the  expected  physiology  range, 

particularly  when  the  //^(r')and  are  small  in  the  background  tissue  regions  or 

are  large  in  the  malignant  cancer  regions.  This  was  exactly  the  problem  that  we  had  with  our 
first  frequency  domain  system  where  the  wavelength  pair  of  (780nm,  830nm)  was  used. 


Fig.3.  The  coefficient 

is  positive 

as  shown  by  circled  line,  however,  the 
coefficient 

is  positive  only  for  wavelength 
<  798  nm  (  see  *  line). 


Fig.  4.  Oxygen  saturation  (Y%  )  vs.  the 
ratio  for  wavelength  pairs  (660nm,  SSOnm, 
circles),  (710  nm,  830nm,  starts)  and 
(780nm,  830nm,  crosses).  Both  pairs  of 
(660nm,  830nm)  and  (710nm,  830)  are 
accq)table. 


Therefore,  to  ensure  equal  sensitivity  of  total  hemoglobin  concentration  estimation  to  two  wavelengths, 
the  wavelength  pair  of  (780  nm,  830  nm)  is  optimal;  while  to  ensure  robust  oxygenation  saturation 
estimation,  the  wavelength  pair  of  (660nm,  830)  is  optimal.  As  a  result,  three  typical  wavelengths  of 
660nm,  780mn,  and  830nm  are  required  for  total  hemoglobin  and  oxygen  saturation. 
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Figure  5  shows  results  obtained  from  a  human  breast  cancer  using  the  optimal  wavelengths  discussed. 
Fig.5.  (a)-(c)  are  ultrasound  images  of  the  breast  cancer  under  chemotherapy  treatment  acquired  before 
treatment,  during  treatment  and  after  six  cycles  of  treatment.  This  example  is  attended  to  demonstrate 
the  use  of  optimal  wavelength  and  the  algorithms  to  calculate  the  total  hemoglobin  concentration  and 
oxygen  saturation,  (d),  (e)  and  (Q  are  corresponding  optical  hemoglobin  concentrations  acquired  at  the 
same  time  as  the  ultrasound  images.  (  g),  (h)  and  (i)  are  oxygen  saturation  maps.  The  shces  in  (d)-(i)  are 
numbered  from  left  to  right  and  top  to  bottom.  Slice  1  is  0.5cm  deep  into  the  breast  and  slice  7  is  3.5  cm 
toward  the  chest  wall.  The  oxygen  saturation  value  is  within  the  physiological  range.  Currently,  we  are 
preparing  manuscript  of  this  work  for  publication. 


Fig.5.  (a)-(c)  ultrasound  images  of  a  human  breast  cancer  under  chemotherapy  treatments,  (d), 
(e)  and  (f)  are  corresponding  optical  hemoglobin  concentrations  acquired  at  the  same  time  as  the 
ultrasound  images(  g),  (h)  and  (i)  are  oxygen  saturation  maps.  The  slices  in  (d)-(i)  are  numbered 
from  left  to  right  and  top  to  bottom.  Slice  1  is  0.5cm  deep  into  the  breast  and  slice  7  is  3.5  cm 
toward  the  chest  wall. 


KEY  RESEARCH  ACCOPLISHMENTS: 

•  Constructed  a  64-channel  receiving  ultrasound  system  for  the  quantification  of  mechanical 
stress  generated  by  the  ultrasound  vibration 

•  Investigated  the  theoretical  basis  of  using  three  wavelengths  to  compute  the  hemoglobin 
concentration  and  oxygenation  from  optical  measurements. 

•  Conducted  in  vivo  study  to  validate  using  three  wavelengths  to  compute  the  hemoglobin 
concentration  and  oxygenation 

REPORTABLE  OUTCOMES: 

•  A  journal  paper  titled  as  “Robust  estimation  of  tumor  hemoglobin  and  oxygenation 
saturation  distributions”  is  in  preparation  based  on  the  work  reported  here. 

CONCLUSIONS: 

I  have  successfully  finished  Task  2  specified  in  the  statement  of  my  proposal.  I  have  built  a  64-chaimel 
receiving  ultrasound  system  to  optimize  the  ultrasound  instrument,  and  have  developed  detailed  analysis 
and  algorithm  on  robust  estimation  of  tumor  Hemoglobin  and  oxygenation  saturation  distributions. 
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